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scattering both in the elastic and inelastic channels and show that the cross section can be expressed 
in terms of the polarized equivalent photon distribution of the proton. We provide the necessary 
kinematical constraints to extract the polarized photon content of the proton using this process 



at HERMES, COMPASS and eRHIC. We also discuss the suppression of the major background 
process coming from virtual Compton scattering. We point out that such an experiment can give 



valuable information on gi(xB,Q 2 ) in the small xb, broad Q 2 region at the future polarized collider 



^ ■ eRHIC and especially in the lower Q 2 , medium xb region in fixed target experiments. 
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I. INTRODUCTION 



QED Compton process (QEDCS) in the scattering Ip — > I'yX, where I is a lepton, has 
a distinctive experimental signature: both the outgoing lepton and photon are detected at 
large polar angles and almost back to back in azimuth, such that their transverse momenta 
almost balance each other, with little or no hadronic activity at the detector . In fact, 
such a reaction in the unpolarized e — p scattering has long been suggested as an excellent 
channel not only to determine the structure function Fo^Xb^Q 2 ) of the proton but also 
to extract the equivalent photon content of the proton (jj, |2j, y|. In a recent Monte Carlo 
analysis of the QED Compton process performed by some members of the HI collaboration at 
HERA j^J, it was found that, although the cross section in the elastic channel was accurately 
described by the equivalent photon approximation (EPA), this was not the case in the 
inelastic channel. In two previous papers 0, 0| we have suggested improved kinematical 
cuts for a more accurate extraction of the unpolarized equivalent photon distribution of 
the proton which furthermore suppress the major background process coming from virtual 
Compton scattering (VCS). In this work we study the QED Compton process in the polarized 
scattering Ip — > /7X (both elastic and inelastic channels), where the initial lepton and proton 
are longitudinally polarized. We show that when the virtuality of the exchanged photon is 
not too large, the cross section can be approximated as a convolution of the longitudinally 
polarized equivalent photon distribution of the proton [?], ls| and the real photoproduction 
cross section. We provide the necessary kinematical constraints to extract the polarized 
photon content of the proton at HERMES, COMPASS and eRHIC (the future polarized 
ep collider planned at BNL). In addition, we show that such an experiment can also access 
the polarized structure function gi(xB, Q 2 ) at HERMES in the low Q 2 region and at eRHIC 
over a wide range of the Bjorken scaling variable xb and Q 2 . gi{xs,Q 2 ) and its first xb 
moment in the low Q 2 region have been studied in fully inclusive measurements at SLAC 
Q, HERMES 0, [n| and JLab 12, Q. The most recent measurements by CLAS Q 



are in the kinematical region Q 2 = 0.15 - 1.64 GeV 2 . The low Q 2 region is of particular 
interest because contributions due to nonperturbative dynamics dominate here and thus 



the transition from soft to hard physics can be studied. In fact the data in 14| clearly 
indicate a dominant contribution from the resonances and at higher Q 2 it is below the 
perturbative QCD evolved scaling value of g%. This in fact illustrates the necessity of further 
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investigation of gi(xB,Q 2 ) in the transition region. In these fixed target experiments, low Q 2 
is associated with low values of Xb, thus the covered kinematical region is smaller compared 
to the unpolarized data. Data on gi{xB,Q 2 ) for small xb and in the scaling region are 
missing due to the absence of polarized colliders so far (with the exception of RHIC, which 
has started operating in the polarized mode for pp collisions only very recently). The small 
xb region is again interesting; it is the region of high parton densities, and measurements 
in this region will provide information about the effects of large (a s ln 2 -^) k resummation 
and DGLAP evolution, and also about the 'soft' to 'hard' scale transition Q,Q,Q. a 
better understanding of gi(xs, Q 2 ) in this region is necessary in order to determine its first 
moment experimentally. The kinematics of QED Compton events is different from the one 
of inclusive deep inelastic scattering due to the radiated photon in the final state and thus 
it provides a novel way to access gi(xs,Q 2 ) in a kinematical region not well covered by 
inclusive measurements (also for F 2 (xb,Q 2 ) ij)- 

The plan of the paper is as follows. In section II and III, we derive the analytic expressions 
of the cross section for the polarized QED Compton process in the elastic and inelastic 
channel, respectively. In section IV, we discuss the background coming from virtual Compton 
scattering (VCS) and also the interference between QEDCS and VCS. The numerical results 
are presented in section V. Summary and conclusions are given in section VI. The analytic 
expressions of the matrix elements are given in Appendices A and B. 



II. ELASTIC QED COMPTON SCATTERING 

We consider QED Compton scattering in the elastic process: 

e(l) + p(P) - e(l')+j(k') +p(P'), (2.1) 

where the incident electron and proton are longitudinally polarized and the 4-momenta of 
the particles are given in brackets. Instead of the electron, one can also consider a muon 
beam (COMPASS); the analytic expressions will be the same. We introduce the invariants 

S=(P + l) 2 , s=(l + k) 2 } t = k 2 . (2.2) 

k = P — P' is the 4-momentum of the virtual photon. The photon in the final state is 
real, k' 2 = 0. We neglect the electron mass everywhere except when it is necessary to avoid 
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divergences in the formulae and take the proton to be massive, P 2 = P' 2 = m 2 . The relevant 
Feynman diagrams for this process are shown in Fig. 1, with X being a proton and Px = P' ■ 
The squared matrix element can be written as 

| AM2 EDCS I 2 = ^H^(P, P')T*(l, k; k'), (2.3) 

H^ U (P,P') and T^(l, k\V , k') being the antisymmetric parts of the hadronic tensor and 
leptonic tensor respectively, which contribute to the polarized cross section. As before 0, 
we use the notation 

/ N.N J3p. 

dPS N {P- P U P N ) = ^f5[P - £ P t j n j^y^o (2-4) 

for the Lorentz invariant iV-particle phase-space element. The cross section can be written 

as 

^oi(S) = 2{s l _ m2) J dPS 2+1 (l + P;l',k',P')\ AM2 EDCS I 2 • (2.5) 
Following the same approach as in we can write this as 

Aa cl (S) = 2{s _ m2) J ^ dPS 2 (l + P-V + k', P') T Xr{P, P')X^{h k) . (2.6) 
Xfi v contains all the informations about the leptonic part of the process and is defined as 

X*(l, k)= J dPS 2 (l + k; I', k')T£,(l, k; l', k'), (2.7) 
T^ u (l, k; I', k') is the antisymmetric part of the leptonic tensor, 

4?> 4 r i 

T^(l,k;l',k') = -—e^kJ{s-t)l a + {u-t)r a }. (2.8) 
su 1 J 

Here e 2 = Ana and we have defined t = (I — I') 2 and u = (/ — k') 2 . 

For polarized scattering, X^ v is antisymmetric in the indices fi, v and can be expressed 
in terms of the Lorentz scalar X 2 : 

X a = -jj^e^k a l Xf(s, t), (2.9) 

with 

Xifrt) = -2X*PZ . (2.10) 



4 



P% is the antisymmetric part of the photon polarization density matrix 



-pliv 



(e'V* - e v e' J *) 



2J\k 2 



(2-11) 



where t a is the spin vector of the photon: 

k 2 1 
t a = N t [k a - —l a )- N t = ^= 



(2-12) 



We define the functions X A (s, t,t) as 



X?(s,t)=2nl dtXf(s,t,t). 



(2-13) 



The integration limits are given by Eq. (2.18) of [5j]. X A (s, t,t) can be obtained using the 
leptonic tensor given above: 

4a 2 f,„ ,2 2ti{u-t) 



x£(s,t,t) 



su(s — t) 



(§ - ty + 



s - t 



-(s + ty 



(2-14) 



The hadronic tensor for polarized scattering is expressed in terms of the proton form 
factors as hi |20|1 , 



H± = -ie 2 e^ pa mk r 



(IV 



nn n q G m {Gm - G E ) k ■ S 

Akj E \j jv/o ct : —r n 



(2.15) 



1 + r m 2 

with t = ^2 and S a = ^(P a — jbla) being the spin vector of the proton which satisfies 
S 2 = — 1 and P ■ S — 0. G E and Gm are the proton electric and magnetic form factors and 
are empirically parametrized as dipoles: 

1 



G E (t) 



[1 -t/(0.71GeV 2 )] 2 
The elastic cross section can then be written as 



, G M (t) = 2.79 G E (t). 



(2.16) 



A 



a 



Ocl 



1 



87r(S-m 2 ) 2 
S-m 2 
s-t 

21^-1 
s-t 



/■'max flf fW . /"27T 

ds —I dt d<f)X A (s,t,t) 

Strain v " ^min ** 

2m 2 s — t 



G 2 M (t) 



t S — m? , 

m 2 s-t \ G M (G M - G E ) 
t S — m 2 



1 + r 



(2-17) 



cf) is the azimuthal angle of the outgoing e — 7 system in the center of momentum frame. 
The limits of integration are the same as in Eqs. (2.18) and (2.24) of ^]. These limits are 
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modified due to the experimental cuts which we impose numerically. In the EPA limit, we 
neglect 1 t | vs. s and m 2 vs. S and get 

Xf{s,t,t) a Xf(s,0,t) = —(---) = -j , 2.18 

S \U S J 71 dt 

where ^ is the differential real photoproduction cross section and 0) = — 4sA<r. The 

elastic cross section then becomes 

a fpa f {l ~Vs^ 2 , f° ,* A , ,dAa(xS,t) 
Aa e i^Aa^ PA =/ tte / rftA 7cl (x) V-^ (2.19) 

where m e is the mass of the electron and A^ e ](x) is the elastic contribution to the polarized 
equivalent photon distribution of the proton jj| 



A 7e i(x) =-^J J (2-x + - 2(1 - x + 



2m 2 x 2 \ 2 / m 2 x 2 \G M {GM - G E 



1 + T 



(2.20) 

with x = -§ and the limits of integration are given by Eq. (2.30) of 0|. 



III. INELASTIC QED COMPTON SCATTERING 

We next consider the corresponding inelastic process 

e(l)+p(P) -> e(J0 +7(*0 +*(P*), (3-1) 

where Px = Px* is the sum over all momenta of the produced hadronic system. We 
take the invariant mass of the produced hadronic system to be W. The Bjorken variable xb 
is defined as 

Q 2 _ Q 2 
%B ~ 2P-{~k) ~ Q 2 + W 2 -m 2 ' (3 ' 2) 

where Q 2 = -k 2 = -t. The cross section for inelastic scattering reads: 

= iw* a -grwrxt, (3.3) 

V / mm e ^ mm W 

where is given by Eq. ()2.9|) and Wj^" is the hadronic tensor: 



= te 2 -^-e^k p 

1 ' K 



1 Q 

9l (k 2 , P ■ k)S a + g 2 (k 2 , P-k)(S a - —P a 



(3.4) 



6 



with S a = —(.Per — jr^a) being the polarization of the proton, satisfying S 2 = — 1 and 
P-S = 0. 

The cross section can be written as 



Aff ine l(5) = ^ IdW 2 ds ' 



a 



A7r(S-m 2 ) 2 Jwl in Jmi J Q l in Q 2 (W 2 + Q 



2 > > ) 2 ' 



X 



m 

S-m 2 W 2 + Q 2 -m 2 2m 2 ( s + Q 2 



2 s + Q 2 + Q 2 + O 2 V S - n, 1 

4m 2 
W 2 + Q 2 - m 



2\ 



9i(x B ,Q 

+ , 92(x B , Q 2 ) Q 2 ), (3.5) 



here X£(s,Q 2 ) = 2tt diX^(s, Q 2 , i) with X 2 A (s,g 2 ,t) given by Eq. (J2~T4]) . The limits 
of the Q 2 , W^ 2 and i integrations are given in Eqs. (3.11), (3.12) and (2.18) of p respectively. 
In the limit of the EPA, as before, we approximate S — m 2 rs S and s + Q 2 ~ s; the cross 
section then becomes 

Aa inel (S) « A^ A = / dx / dt A 7inel (x,xS) (3.6) 

Jx min Jm\-s at 

where again x = s/S and A7 ine i(x, xS") is the inelastic contribution to the polarized equiva- 
lent photon distribution of the proton: 

A 7mcl (x,^) = -^ -y^ -Q^{ 2 -y--Q^)H? Q h (3 - 7) 

where we have taken the scale to be s. Here we have neglected the contribution from 
92{xb,Q 2 )- Expressing gi(xB,Q 2 ) in terms of the polarized quark and antiquark distribu- 
tions, one can confirm that the above expression reduces to that given in [l,^. However, in 
this case, one chooses the minimal (but not compelling) boundary condition A 7 (x, Qq) = 
at a scale Ql = 0.26 GeV 2 . Eq. ()3.7j) is free from this particular boundary condition. The 
limits of the Q 2 integration can be approximated similar to the unpolarized case (see Eq. 
(3.16) of U). 



IV. BACKGROUND FROM VIRTUAL COMPTON SCATTERING 

The cross section of the process in Eq. ()3.1j) (also the elastic channel) receives contribution 
from the virtual Compton Scattering (VCS), when the photon is emitted from the proton 
side (see Fig. 2) as well as the interference between the QED Compton scattering (QEDCS) 
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and VCS. The cross section for the elastic process is given by 



A(J el(S) = ^ / ds dt I dt / d<j)- -| AM e i \, (4.1) 

Sn(b — m 1 ) 1 J ml Jt min Jt min Jo [s — t) 

where 

| AM cl | 2 = | AM2 EDCS | 2 + | AM^ CS | 2 - 2 AM^ EDCS AMl cs * (4.2) 

is the matrix element squared of the subprocess. The limits of integrations are the same as in 
Eq. (j2.17|) . The interference term will have opposite sign if we consider a positron instead of 
an electron. The cross section of the VCS process is expressed in terms of generalized parton 
distributions and one needs a realistic model for a quantitative estimate of this background 
j21|. Here, in order to find the constraints to suppress the VCS, we make a simplified 
assumption: we take the proton to be a massive pointlike particle with an effective 7*p 
vertex, —i^Fi(t). The explicit expressions for the matrix elements are given in Appendix 
A. 

Particularly interesting for our purpose of extracting the polarized photon distribution 
of the proton is the inelastic channel. Here we use a unified parton model (similar to our 
analysis in j^]) to estimate the VCS and QEDCS rates. The cross section within the parton 
model is given by 

dAa- ine \ ^ . , dAa q 



TAq(x B ,Q 2 ) " s , (4.3) 

dxB dQ 2 ds dt dip q ds dQ 2 dt dcf) 

where Aq(xB,Q 2 ) are the polarized quark and antiquark distributions of the initial proton, 
q — u, d, s, u, d, s and dAa q is the differential cross section of the subprocess 

eil) + qip) - e(0 + 7 (*0 + q(p)- (4-4) 

Here q is a longitudinally polarized quark in a longitudinally polarized proton and q is a 
quark in the final state. The integrated cross section becomes 

A(T inel (S) = — — / dW 2 / ds / / dt / d6- — 

meU ; %n{S - m 2 ) 2 ^ Jw^ Jml Jq^ Q 2 Ji min Jo (s + Q 2 ) 

x | AM incl | 2 Aq(x B ,Q 2 ), (4.5) 

where 

| AM incl | 2 = | AM?J DCS | 2 + | AM^ S | 2 - 2 VteAM?J DCS AM™** . (4.6) 



The limits of integrations are the same as in Eq. (17) of The explicit expression of 
the matrix element is given in the appendix B. It is useful to define the auxiliary invariants 
S = (p' + k') 2 and U = (p' — k) 2 , which can be written in terms of measurable quantities, 

S= ^ Xl - XB \ (j = i-S + Q 2 , (4.7) 

with xi = 2 p\i-i>) ■ m addition to the leptonic variable x\ we define x 7 = -p^, which represents 
the fraction of the longitudinal momentum of the proton carried by the virtual photon jsj]. 
In the limit of the EPA, both x\ and x 7 are the same and become equal to x. 

V. NUMERICAL RESULTS 

In this section, we present our numerical results. The cuts used for HERMES, COMPASS 
and eRHIC kinematics are given in Table I. The constraints on the energies and polar angles 
of the detected particles reduce the background contributions coming from the radiative 
emissions (when the final state photon is emitted along the incident or the final lepton line), 



i. 



2J. QED Compton 



because they prevent the lepton propagators to become too small 
events are singled out at HERA by imposing a maximum limit on the acoplanarity angle 
<ft a which is defined as <pA = [7T — | 7 — <p e | |, where 7 and (f) e are the _azimuthal angles 



of the outgoing photon and electron, respectively. We have observed p, |6( that instead of 
this limit on <pA, the constraint s > Q 2 , which is applicable experimentally, is more efficient 
in extracting the equivalent photon distribution from the 'exact' result. Here we use this 
constraint. 

i—i 

The unpolarized cross sections have been calculated using the formulae in |5[ but for 
HERMES, COMPASS and eRHIC kinematics, respectively. In the numerical estimate of 
the unpolarized cross section, we have used ALLM97 parametrization j^j of the structure 
function ^(ib, Q 2 ) as before, which is obtained by fitting DIS data of HERA and fixed target 
experiments together with the total pp and 'jp cross sections measured and is expected to 
hold over the entire range of xb and CP. We have taken Fl(xb, Q 2 ) to be zero, assuming the 
Callan-Gross relation, similar to % |22|. In the polarized cross section, we have neglected 

n 

the contribution from g2(xB,Q 2 ) and used the parametrization of 24] for gi(xB,Q 2 )- In 
this parametrization, gi{xB,Q 2 ) is described in the low-Q 2 region by the GVMD model 
together with the Drell-Hearn-Gerasimov-Hosoda-Yamamoto sum rule and the asymptotic 
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part of gi(xB, Q 2 ) is expressed in terms of NLO GRSVOO [25[ parton distributions (standard 
scenario) in terms of a suitably defined scaling variable x = q 2+ q2 +W 2_ M 2 with Qq — 1.2 
GeV 2 . The scale Q 2 is changed to Q 2 + Q 2 ,, so as to extrapolate to low-Q 2 region. It is to 
be noted that for QED Compton scattering, the effects of Fl(xb,Q 2 ) and #2(^5, Q 2 ) have 
to be taken into account in a more accurate study as their effect may become non-negligible 
in the low-Q 2 region. However, this is beyond the scope of the present work. 

Before discussing the results for specific experiments, it is interesting to investigate some 
general properties of the total cross section. Figs. 3 (a) and (b) show the total QEDCS cross 
section, polarized and unpolarized, respectively, as a function of the incident lepton energy 
Ei. We have imposed the constraints in the second column of Table I on the enegies and 
angles of the outgoing particles, as well as those on s. Both polarized and unpolarized cross 
sections increase sharply with Ei, reach a peak at around E\ = 20 GeV and the start to 
decrease. The cross section in the inelastic channel is also shown, which have similar trends 
except that the peak in the polarized case is broader. 

A. HERMES 

Figs. 4(a) and (b) show the total (elastic+inelastic) polarized and unpolarized QED 
Compton scattering cross sections, respectively, in bins of x 7 for HERMES kinematics, 
subject to the cuts of Table I. We have taken the incident electron energy E e = 27.5 GeV. 
We also show the cross section calculated in the EPA. The same in the inelastic channel is also 
shown. The cross section, integrated over x 7 , agrees with the EPA within 7.1% (unpolarized) 
and 4.8% (polarized). From the figures it is also clear that the agreement in the inelastic 
channel is much better than for HERA kinematics j^lO] (2.5% in the polarized case). This is 
because at HERMES Q 2 can never become too large (maximum 13.7 GeV 2 ), subject to our 
kinematical cuts, which is expected in a fixed target experiment. The agreement is not so 
good without the constraint s > 1 GeV 2 . Fig. 4(c) shows the asymmetry, which is defined 
as 

A LL = 5.1 

where + and — denote the helicities of the incoming electron and proton. They are calculated 
with the same set of constraints. The asymmetry is quite sizable at HERMES and increases 



10 



in higher x 7 bins. The asymmetry in the EPA is also shown. It is interesting to note 
that the discrepancy in the cross sections with the EPA approximated estimate, actually 
gets canceled in the asymmetry, as a result it shows an excellent agreement with the EPA, 
except in the last bin. We have also shown the expected statistical error in the bins, which 
have been calculated using the formula, valid when the asymmetry is not too large: 

SA LL « 1 ; (5.2) 

where V e and V v are the polarizations of the incident lepton and proton, respectively, C is 
the integrated luminosity and aun is the unpolarized cross section in the corresponding x 7 
bin. We have taken V e = V p = 0.7 and £ = lfb^ 1 for HERMES. The expected statistical 
error increases in higher x 7 bins, because the number of events become smaller. However 
the asymmetry seems to be measurable at HERMES. 

The background from virtual Compton scattering is reduced at HERA by the experi- 
mental condition of no observable hadronic activity at the detectors. Basically the electron 
and photon are detected in the backward detectors and the hadronic system in the forward 
detectors. In our previous work, we have observed that for unpolarized scattering at HERA, 
such a constraint is insufficient to remove the VCS contribution for higher x T We have 
proposed a new constraint S > s, where S and s can be measured experimentally, to be im- 
posed on the cross section. Here, we investigate the effect of this constraint on the polarized 
cross section. To estimate the inelastic contribution coming from VCS, we use Eq. (|4.5j) . 
together with an effective model for the parton distribution of the proton. The effective 
parton distribution is of the form 

Aq(x B ,Q 2 ) = Aq(x,Q 2 + Ql), (5.3) 

Aq(x B ,Q 2 ) being the NLO GRSV00 (standard scenario) distribution function [25]. In the 
relevant kinematical region, Q 2 can be very small and may become close to zero, where 
the parton picture is not applicable. The parameter Qf, = 2.3 GeV 2 prevents the scale 
of the parton distribution to become too small, x is a suitably defined scaling variable, 

= _ SB(Q 2 +<ffl 

X Q 2 +x B Ql ■ |-| 

To estimate the unpolarized background effect, we use the same expressions as in p with 
an effective parton distribution given in Eq. (22) of Fig. 4(d) shows the polarized cross 
section in the inelastic channel at HERMES, subject to the constraints of Table I, in bins of 
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s — S calculated in the 'effective' parton model. The VCS and the interference contributions 
are also shown. QEDCS cross sections using the Badelek et. al parametrization of gi(xs, Q 2 ) 
are also plotted. In fact, the cross section in the 'effective' parton model lies close to this. 
Within theparton model, VCS is suppressed when s < S, similar to the unpolarized case 
at HERA |6j|. Unlike HERA, the interference between QEDCS and VCS is not negligible at 
HERMES, although smaller than QEDCS in the relevant region. Since the interference term 
changes sign when a positron beam is used instead of the electron beam, a combination of 
electron and positron scattering data can eliminate this contribution. In order to estimate 
the VCS in the elastic channel, one needs a suitable model for the polarized generalized 
parton distributions. However, in the simplified approximation of a pointlike proton with 
an effective vertex as described in section IV, the elastic VCS as well as the interference 
contribution is much suppressed at HERMES. Similar observations hold for unpolarized 
scattering. 

Fig. 5 shows the asymmetries in the inelastic channel in bins of xb- In addition to the 
cuts mentioned above and shown in table I, we have also chosen S — s > 2 GeV 2 to suppress 
the background. The asymmetry is small but sizable and could be a tool to access gi(xs, Q 2 ) 
at HERMES. In fact, QED Compton events can be observed at HERMES in the kinematical 
region xb = 0.02 — 0.7 and Q 2 = 0.007 — 7 GeV 2 (small Q 2 , medium xb)- However, from 
the figure it is seen that the asymmetry is very small for xb below 0.1. We have also shown 
the expected statistical error in each bin. The average Q 2 value in GeV 2 for the polarized 
cross section for each bin is shown, which has been calculated using the formula 

Jbin dA(J 

B. COMPASS 

Figs. 6 (a) and (b) show the cross sections in bins of x^ for the polarized and unpolarized 
QEDCS for the kinematics of COMPASS. We take the energy of the incident muon beam 
to be 160 GeV, the target is a proton. The final muon and the photon are detected in 
the polar angle region 0.04 < 9 IX ,9 1 < 0.18. The cross sections in bins, subject to the 
kinematical constraints shown in Table I, are much smaller than at HERMES, because they 
start to decrease with the increase of the incident lepton energy Ei as E\ becomes greater 
than about 20 GeV, as observed in Fig. 3. As before, the cuts remove the initial and final 

12 



state radiative events. The x 1 integrated cross section agrees with the EPA within 14.2% 
(unpolarized) and 15.5% (polarized). The agreement thus is not as good as at HERMES. 
From the figures it is seen that the coss section in the EPA actually lies below the 'exact' 
one, both for polarized and unpolarized cases. This discrepancy is due to the fact that the 
EPA is expected to be a good approximation when the virtuality of the exchanged photon is 
small. At COMPASS, with our kinematical cuts, Q 2 can not reach a value below 0.07 GeV 2 
and can be as large as 144 GeV 2 , whereas for HERMES smaller values of Q 2 are accessible 
(see the previous subsection). Fig. 6(c) shows the asymmetry in bins of x 7 , the asymmetry 
in the inelastic channel is also shown. The asymmetry is of the same order of magnitude 
as in HERMES and is in good agreement with the EPA. We have also shown the expected 
statistical error in each bin, calculated using Eq. ([5.2)1 . We have taken V e = V p = 0.7 
and C = lfb" 1 for COMPASS. The statistical error is large in higher x 7 bins. Fig. 6(d) 
shows the polarized QEDCS, VCS and interference contributions (inelastic) calculated in 
the 'effective' parton model, in bins of s — S. As in HERMES, VCS is suppressed for s < S. 
The interference term is not suppressed but using /i + and jjT beams this can be eliminated. 
We have also plotted the QEDCS cross section using Badelek et al. parametrization of 
9i{xb,Q 2 )- The VCS and the interference contributions (elastic) are much suppressed in 
the pointlike approximation of the proton with the effective vertex. Fig. 7 shows the 
asymmetry at COMPASS in the inelastic channel plotted in bins of xb with the same set of 
constraints and the additional cut S — s > 2 GeV 2 . The asymmetry is sizable and can give 
access to gi(xB,Q 2 ), the kinematically allowed range is 0.07 < xb ■ We have also shown 
the expected statistical errors in the bins and the average Q 2 in each bin. Comparing Fig. 
5 and 7 one can see that there is no overlap in the kinematical region covered at HERMES 
and COMPASS. Higher values of Q 2 are probed at COMPASS in the same xb range as 
compared to HERMES. 



C. eRHIC 



The cross sections for eRHIC kinematics, both polarized and unpolarized, are shown in 
Fig. 8(a) and (b), respectively, in bins of x T We have taken the incident electron energy 
E e = 10 GeV and the incident proton energy E p = 250 GeV. The cross section in the EPA 
is also shown. The kinematic constraints are given in Table I. The polar angle acceptance of 
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the detectors at eRHIC is not known. We have taken the range of 8 e , # 7 to be the same as 
at HERA. We have checked that the constraints on the energies and the polar angles of the 
outgoing electron and photon are sufficient to prevent the electron propagators to become too 
small and thus reduce the radiative contributions. The unpolarized total (elastic+inelastic) 
cross section, integrated over x 1 agrees with the EPA within 1.6%. The agreement in the 
inelastic channel is about 6.3%. The polarized total cross section agrees with the EPA within 
9.8%. The EPA in this case lies below the 'exact' one in all the bins. The agreement in 
the inelastic channel is about 19.6%. More restrictive constraints instead of s > Q 2 , like 
s > 10 Q 2 , makes the agreement better, about 1.2% in the polarized case and 1.9% in the 
unpolarized case. Fig. 8(c) shows the asymmetry for eRHIC, in bins of x 1 . The discrepancy 
in the cross section cancels in the asymmetry, as a result good agreement with the EPA 
is observed in all bins except the last one at higher x T The asymmetry in the inelastic 
channel is also shown. We have plotted the expected statistical error in the bins using Eq. 
(|5.2J) . For eRHIC, we have taken V e = V p = 0.7 and £ = lfb^ 1 . The expected statistical 
error increases in higher x 1 bins. The asymmetry is very small for small x 7 but becomes 
sizable as x 1 increases. Fig. 8(d) shows the polarized cross section in the inelastic channel, 
in bins of s — S, in the 'effective' parton model for eRHIC. The VCS is suppressed in all 
bins, especially for s < S. The interference contribution is negligible, similar to HERA. The 
'effective' parton model QEDCS cross section is also compared with the more exact one, 
using Badelek et. al. parametrization for gi(xB,Q 2 )- Similar effects are observed in the 
unpolarized case. In the pointlike approximation of the proton with the effective vertex, as 
before, the elastic VCS as well as the interference contributions are very much suppressed. 
Fig. 9 shows the asymmetry in bins of xb in the inelastic channel, which may be relevant 
for the determination of gi(xB,Q 2 ) using QEDCS at eRHIC. The asymmetry is small but 
sizable, however the error bars are large and therefore good statistics is needed, xb can be 
as low as 0.002. A wide range of Q 2 can be accessed at eRHIC starting from 0.008 to 2000 
GeV 2 ; the average Q 2 value in the bins ranges from 2.4 to 315 GeV 2 . Fig. 10 shows the 
total asymmetry in Q 2 bins for eRHIC. The asymmetry in this case is bigger in each bin 
and the error bars are smaller than the xb bins except the last bin for high Q 2 where the 
number of events are smaller. 
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VI. SUMMARY AND CONCLUSIONS 



To summarize, in this paper we have analyzed the QED Compton process in polarized 
Ip scattering, both in the elastic and inelastic channel. This process has a distinctive ex- 
perimental signature and we showed that the cross section can be expressed in terms of 
the equivalent photon distribution of the polarized proton, convoluted with the real pho- 
toproduction cross section. The EPA is a useful tool to estimate high energy scattering 
cross sections; however the accuracy of this approximation and the kinematical region of 
its validity have to be checked by experiment. In this work we provided the necessary 
kinematical constraints for the extraction of the polarized photon content of the proton by 
measuring the QED Compton process at HERMES, COMPASS and eRHIC. We showed 
that the cross section and, in particular, the asymmetries are quite accurately described by 
the EPA. We also discussed the possibility of suppressing the major background process, 
namely the virtual Compton scattering. We pointed out that such an experiment can give 
access to the spin structure function gi(xB,Q 2 ) in the region of low Q 2 and medium xb in 
fixed target experiments and over a broad range of xb, Q 2 at the future polarized ep collider, 
eRHIC. Because of the different kinematics compared to the fully inclusive processes, the 
QED Compton process can provide information on gi(xs,Q 2 ) in a range not well-covered 
by inclusive measurements and thus is a valuable tool to have a complete understanding of 
the spin structure function. 
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APPENDIX A: MATRIX ELEMENT FOR THE ELASTIC BACKGROUND VIR- 
TUAL COMPTON SCATTERING PROCESS 



The explicit expressions of the matrix elements in section IV are given below: 



AM QEDCS 



tsii 



~ A + T¥ B . 



(Al) 



with 



AM 



VCS I 2 



cl 



tU' S> 



2m2 n 

A + C 

5' S'U' . 



A = 2t 2 + (s - 2S' -u)(s + u)- 2t(s - 2S' - U') - 2uU' , 



(A2) 



(A3) 



B = -2t 6 + s- su 2 + 2r(2s + u)- t(3s 2 + it 2 ), 



(A4) 



C = (s + u) 2 [-2S' 2 + su- 2Su - u 2 + 2m 2 (s + u - t) - t(S' + «)] 
-(s + u) [2t 2 - 3ts + s 2 + s(S' - 2u) + 3u(S' + u)] U' 
- [2t 2 + s 2 -su + 2u 2 - t(3s + it)} U' 2 . 



2?HeAM2 EDCS AM^ cs * 



4 



tsutU' S' 



D 



2m 2 E 



F 1 (t)F 1 (t), 



where D and E read: 



D = [2t 2 + (s - 2S' -u)(s + u)- 2t(s - 2S' - U') - 2uU'\ 
x {(s + u) [tii + S'(s + u)] + [t(s -u) + s(s + u)]U'}, 



(A5) 
(A6) 



(A7) 



E = [s(s - t) 2 (s - 2t) + (2t 3 - t 2 s - s 3 )u + {-2t 2 - 3ts + s 2 )u 2 + (t + 3s)u 3 ] U' 
-(s + u){-2t 3 u - (s + u)[s 2 (S' - 2u) + S'u 2 + 2su{S' + u)} 
-t[-7sS'u + u 2 {-2S' + u) + s 2 (-3S' + 5m)] + t 2 [2u{u - S') + s{5u - 25")]}. 

(A8) 

We have introduced the invariants U = (P — k') 2 , u = (l — k 1 ) 2 and 5 = — (s+u+U' — m 2 ) 
and used the notations S' = S — m 2 , U' = U — m 2 , 5" = 5 — m 2 for compactness. 
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APPENDIX B: MATRIX ELEMENT FOR THE INELASTIC VIRTUAL COMP- 
TON SCATTERING BACKGROUND PROCESS 

For the corresponding inelastic channel in section IV the explicit matrix elements are 
given by: 

\AM?J D ° S I' = 4e^. (Bl) 

\AM™*f = -4e',j§fi, (B2) 
with S = — (s + u + x B U'), U = x B U' and 

F = s 2 - u 2 + 2Q 4 + 2Q 2 s - 2x B [sS' + u(S' + U') + Q 2 (2S' + [/')]■ ( B 3) 
Here e q is the charge of the parton in units of the proton charge. Finally 

2MeAM?J DCS AM™s* = _ 4e 3 (B4) 

mcl q Q*SUtSU 

with 

G = 2Q 4 + s 2 - 2sS'x B - u[u + 2(S' + U')x B ] + 2Q 2 [s - (2S' + U')x B ] (B5) 
H = Q 2 [s(u — U) + u(u + U)}. -x B (s + u)[S'u + s(S' + U')] (B6) 
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y 




+ Crossed. 



Fig. 1: Fcynman diagrams for the QED Compton process 
(QEDCS). X = p (and Px = P') corresponds to elastic scat- 
tering. 




Fig. 2: As in Fig. 1 but for the virtual Compton scattering (VCS) 
background process. 
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Fig. 3: QEDCS cross section vs. energy of the incident lepton; (a) polarized, (b) unpolarized. 
The continuous line is the total cross section and the dashed line is the cross section in the 
inelastic channel. The cuts imposed are given in the central column of table I. We have used 
the ALLM paramtrization of F2 [3] and the Badelek et. al parametrization of 31 | . 
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Fig. 4: Cross section for QED Compton scattering (QEDCS) at HERMES in bins of x 1 (a) 
)olarized, (b) unpolarized, (c) the asymmetry; for the polarized cross section Badelek et al. 

parametrization of g\ (BKZ) and for the unpolarized cross section ALLM parametriza- 
tion of F2 have been used; (d) polarized inelastic cross section for QEDCS (long dashed), 
VCS (dashed-dotted) and the interference (dashed) at HERMES in the effective parton 
model. The bins are in s — S, expressed in GeV 2 . The continuous line is the QEDCS cross 
section using the BKZ parametrization of giixs, Q 2 )- The constraints imposed are given in 
Table I. 
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Fig. 5: Asymmetry in the inelastic channel in bins of xb at HERMES. We have used Badelek 
et al. 24] parametrization of g\. The constraints imposed are as in Table I (except s > Q 2 ), 
together with S — s > 2 GeV 2 . The average Q 2 (in GeV 2 ) of each bin is also shown. 
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Fig. 6: (a), (b), (c) and (d) are the same as in Fig. 4 but for COMPASS. The constraints 
imposed are given in Table I. 
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Fig. 7: Asymmetry in the inelastic channel in bins of at COMPASS. We have used 
Badelek et al. 24i] parametrization of g\ . The constraints imposed are as in Table I (except 
s > Q 2 ), together with S — s > 2 GeV 2 . The average Q 2 (in GeV 2 ) of each bin is also 
shown. 
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Fig. 8: (a), (b), (c) and (d) are the same as in Fig. 4 but for eRHIC. The constraints 
imposed are given in Table I. 
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Fig. 9: Asymmetry in the inelastic channel in bins of xb at eRHIC. We have used Badelek 
et al. 24 parametrization of g\. The constraints imposed are as in Table I (except s > Q 2 ), 
together with S > s. The average Q 2 (in GeV 2 ) of each bin is also shown. 
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Fig. 10: Asymmetry in bins of Q 2 (GeV 2 ) at eRHIC. We have used Badelek et al. 
parametrization of g±. The constraints imposed are the same as in Fig. 9. 
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HERMES 


COMPASS 


eRHIC 






.Bp = 250 GeV 


E e = 27.5 GeV 


Ep = 160 GeV 


E e = 10 GeV 


0.04 < e e , 6> 7 < 0.2 


0.04 < Op, 6» 7 < 0.18 


0.06 < e e , e 1 < vr - o.o6 


E' e , E\ > 4 GeV 


fy, K f > 4 GeV 


E' e , E^ > 4 GeV 


s > 1 GeV 2 


s > 1 GeV 2 


1 > 1 GeV 2 


s > Q 2 


s > Q 2 


s > Q 2 



TABLE I: Energies, angular acceptance and kinematical cuts for various experiments. 
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